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有 理 三 次 Hermite 插值 样 条 及 其 逼近 性 质 * 
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摘 ”要 : 提高 插值 曲线 曲面 的 逼近 性 是 计算 辅助 几何 设计 中 的 一 个 重要 问题 . 本 文 构建 了 一 种 带 单 参数 的 
分 段 有 理 三 次 Hermite 插值 样 条 . 讨论 了 该 样 条 的 逼近 性 ， 给 出 了 一 种 提高 插值 曲线 曲面 逼近 性 
的 方法 ， 并 且 给 出 数值 例子 . 结果 表明 ， 对 于 给 定 的 插值 条 件 ， 通 过 选择 合适 的 参数 ， 依 本 文 方 
法 所 生成 的 插值 曲线 曲面 在 逼近 效果 上 好 于 标准 三 次 Hermite 插值 曲线 曲面 . 
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插值 法 是 计算 机 辅助 几何 设计 中 用 于 曲线 曲面 造型 的 一 种 重要 工具 . 标准 的 分 段 三 
次 Hermite 插 值 是 其 中 的 一 个 重要 的 方法 . 但 是 ， 在 插值 条 件 确定 的 情况 下 ， 分 段 三 
次 Hermite 插 值 曲 线 的 形状 完全 确定 的 ， 揪 值 曲线 的 逼近 效果 也 随 之 确定 ， 这 种 插值 被 称 
之 为 确定 性 插值 . 如 何在 插值 条 件 确 定 的 情况 下 提高 逼近 效果 ， 是 计算 机 辅助 几何 设计 中 
的 重要 研究 课题 . 近 些 年 来 ， 不 少 作 者 讨论 了 有 理 三 次 插值 样 条 的 逼近 性 彤 ]， 这 些 插值 样 
条 都 具有 标准 的 三 次 Hermite 插 值 相似 的 性 质 . 上 述 文 献 在 讨论 插值 样 条 的 逼近 性 时 ， 应 
用 Peano-Kernel 定理 ， 证 明了 插值 误差 系数 是 有 界 的 ， 从 而 ， 在 理论 上 ， 当 插值 区 间 hi; 趋 于 
无 穷 小 时 ， 播 值 曲 线 逼 近 于 被 插 函 数 ， 但 实际 应 用 过 程 中 ， 插 值 区 间 hh; 不 可 能 趋 于 无 穷 小 ， 
否则 会 导致 计算 量 的 增加 .至 于 在 相同 的 插值 条 件 下 ， 这 些 样 条 在 逼近 效果 上 能 否 比 标准 三 
次 Hermite 插值 样 条 好 ， 这 些 文献 均 没 有 作 比 较 . 另外 ， 这 些 有 理 形 式 的 插值 样 条 含有 多 个 参 
数 ， 使 得 对 它们 的 有 关 逼 近 性 讨论 变 得 复杂 . 

本 文 提出 一 种 分 段 有 理 三 次 Hermite 插 值 曲线 ， 每 段 只 带 一 个 参数 ， 在 讨论 允 近 性 及 计算 
参数 取 值 时 都 比较 简单 ， 对 于 被 插 函 数 ， 在 插值 区 间 不 趋 于 无 穷 小 的 情况 下 ， 若 参数 取 值 合 
适 ， 生 成 的 插值 曲线 在 允 近 效果 上 好 于 标准 的 三 次 Hermite 插值 曲线 . 数值 例子 表明 ， 这 种 参 
数 取 值 方法 是 有 效 的 . 


2 ”有理 三 次 Hermite AKMX MÉI Ferguson 曲线 
我 们 知道 ， 对 于 给 定 的 节点 a = £o < T1 <- < Irn = b, Whi = Ti+ — Ti, t= a] 1E 
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区 间 [zi, zi41] 上 的 标准 三 次 Hermite 基 函 数 具 有 如 下 形式 
Ft) =1- 3+2, Fi(t)-30?-—20, G;(t)2t—20?-0, Gut) =- +t, (1) 
它们 满足 
Fi(0) = Fa (1)=1, F:(1) = Fi+ı(0) = 0, 
F; (0) = F; (1) = Fiy (0) = Fi} (1) = 0, 
Gi(0) = Gi(l) = Gi} (0) = Gia (1) = 0, 
Gi(0) = Gia (1) =1, Gi(1)= Gi (0) = 0, 
BE FU) + F(t) 2 1, Gi(t) = -GaaQ - t). 
基于 这 组 基 函 数 的 三 次 Hermite 插 值 曲线 具有 C1 连续 性 ， 但 插值 曲线 的 形状 是 固定 不 变 
的 . 带 有 参数 的 有 理 形式 的 插值 样 条 可 以 改变 插值 曲线 的 形状 上 1. 
下 面 构建 具有 三 次 Hermite 揪 值 曲线 性 质 的 有 理 形式 的 插值 样 条 曲线 . 先 构建 有 理 形式 的 


PT E 
定义 1 对 任意 实数 和 ; > -2, 0<t<1, 称 


1+Ai)t? Ait? 
, RFiri(t) = 1F(A;-2)0t-(2—À,)t2 





RGi(t) = xc eoe RG (0 = sen tage n 
为 带 有 参数 X; 的 有 理 三 次 Hermite 基 函数 . 
经 简单 计算 知 有 
RF(0) = RFin(1)-71, RF;(1) = RFi(0) = 0, 
RF;(0) = RF;(1) = RF;,1(0) = RF;,(1) — 0, (4) 


)- 
RGi(0) = RG;(1) = RG;41(0) = RGi+1(1) = 0, 
RGi(0) = RG;,,(1)-1, RG;(1) = RGj,(0)- 0, 


RFE(t) + RF (t) = 1, RG(( = -RG (1 - t). 

由 上 面 性 质 可 知 ， 上 面 的 基 函 数 具 有 与 标准 三 次 Hermite 基 函数 相同 的 特性 .特别 地 ， 
当 A = 2 时 ， 有 理 三 次 Hermite 基 函数 退化 为 三 次 Hermite 基 函数 . 因此 ， 可 用 来 作 两 点 
的 Hemite 插值 ， 对 应 的 带 有 参数 的 有 理 三 次 Ferguson 曲线 定义 为 


RHi(t) = RFi(t)pi + RFiti(t)piti + RGi(t)p; + RGitri(t)pit = 0, (5) 


其 中 pi, poi 和 pi, piyi 为 两 个 插值 端点 和 端点 切 矢 . 在 端点 及 端点 切 矢 不 变 的 情况 下 ， 利 用 
参数 X; 的 不 同 取 值 ， 可 以 得 到 不 同形 状 的 Ferguson 曲线 . 图 1 为 有 理 三 次 Ferguson 曲线 ， 其 
中 虚线 为 标准 的 三 次 Ferguson H. 

从 图 象 上 看 ， 参 数 有 共有 明显 的 几何 意义 ， 即 当 参 数 取 值 大 于 或 小 于 2 时 ， 对 应 的 图 象 
在 三 次 Ferguson 曲线 的 下 面 或 上 面 . 这 样 ， 有 理 三 次 Hermite 插 值 曲线 就 能 比 标准 的 三 
次 Hermite 插 值 曲线 更 能 逼近 或 偏离 被 插曲 线 . 因而 有 理 三 次 Hermite 插 值 曲线 有 更 好 的 
实用 性 . 
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图 1: 有 理 三 次 Ferguson 曲线 


3 有理 三 次 Hermite 插值 曲线 


定义 2 给 定数 据 (zi,yi,di), i =0,1,…,n， 此 处 , a = zo < tı <- < zn 二 6 是 分 
划 点 ， 纪 和 ai 为 分 划 点 zx; 处 的 函数 值 及 一 阶 导 数值 . Whi = tipi 一 Xi, t= exa, Hf 
数 入 ; -2. 称 


RHi(x)|iz, s, i) = y RF;(t) RUE yii RF (t) 
+dihiRGi(t) cT di yh RGiQ(t), i=0,1,...,n, (6) 


为 [a, 如 上 的 分 段 有 理 三 次 Hermite 插 值 样 条 . 其 中 RE,(£), RFE;a(t), RG(t) 及 RG t) AA 
理 三 次 Hermite 基 函数 . 

显然 ， 对 给 定 的 数据 (zi; Vi; di), i= 0, 1, t4, RHi(z) 满足 

RHi(zi) = yi, RH;(a) = di, $20,1,-- ,n. 

容易 验证 ， 当 入 ; = 2 时 ， 即 为 标准 的 分 段 三 次 Hermite 插值 样 条 .因而 有 理 三 次 Hermite 插 
值 样 条 是 三 次 Hermite 插值 样 条 的 一 种 推广 . 

假如 我 们 选择 适当 的 参数 Xi, RH (o) 还 能 够 变 成 C2 连续 的 插值 样 条 . 事实 上 , 令 RH" (x4) 
= RH"(zi—-), i 二 1,2,… ,n 一 1， 可 以 得 到 如 下 连续 性 方程 


h2 [((yi-i — yi) (1 4-3) + hi-1(di-1 + NM-1di)] 
= hal -4)-cA)th(u Xd), i212, ,n-1, (7) 


称 上 式 为 C2? 连续 性 约束 条 件 . 
由 此 ， 可 由 [zo,zij 上 的 Xo 通 过 式 (7) 确定 Ai， 再 由 Xi， 通 过 式 (7) AEA KERE, Z 
段 构造 出 [a, b] 上 的 C? 连续 的 有 理 三 次 Hermite 插值 曲线 . 


4 有理 三 次 Hermite 插值 曲线 /曲面 的 逼近 性 


设 被 插 函 数 f(z) 在 插值 区 间 上 是 C? 连续 可 微 的 ， 对 有 理 三 次 Hermite 插值 曲线 ， 有 如 下 的 
误差 估计 定理 . 
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定理 1 设 f(x) € [a,4], a= zo <zi<…<zn,=b 是 区 间 [a,b] 上 的 一 个 分 划 ，RH(x) 为 
有 理 三 次 Hermite 插值 样 条 ， 对 任意 给 定 的 参数 和 i， 23g € [ri zii] FP, A 
[RIFI = fO - RE G)|| < ssh ion. 
证 明 因为 这 种 插值 具有 局 部 性 质 ， 只 需 考 虑 在 [, 引 的 子 区 间 上 的 误差 即 可 .容易 验证 : 
该 种 插值 对 一 次 多 项 式 精 确 成 立 ， 于 是 当 z e [zi, zifi] 时 ， 利 用 Peano-Kernel 定理 回 ， 有 


RU) - He) - RHG) = [^ FOORE- ràn, 


(x = T) = e = t) [3(zi+ == hi)] ey P(zia as hi), PiLT T, 


—t2(1 -—— t) [3(zi.1 -— hi)] == t (Tipi R hi), r«T« Ti 十 1， 


Re[(z 一 r)+] = | 


p(T), Zi <T <T, 
g(T) E<TX< Tipi 


首先 考察 g(7) 在 [z,ziti] 上 的 性 质 ， 容易 得 到 g(xz) = 2? (t-1)(zii1-2) < 0， 且 由 gq(zi#1) 
= P (ziy — 1) > 0 可 知 g(7) 在 (zx,zi+1) 有 一 个 根 ， 可 求 得 该 根 为 


接 下 来 ， 考 察 p(7) 在 (zi, x) 上 的 性 质 ， 因为 p(7) = (z 一 +) 十 gq(7)， 于 是 p(x) = qlr) < 0， 
H.p(z;) = (x — 2i) + q(zi) = t(1— t)?hi >0. 
同样 可 求 得 p(7) YE (x, zi+1) 的 根 为 





, hi(1 +t) 
SE RE Rar 
于 是 就 有 
IE) -REGI < WN f pare f woar] [7 ame 
(=t)? 
Ur 
四 AC(1— t? 
“的 一 n 20 — 25 
易 求 得 
EW = Te 


证 毕 
HIW, SEKE h — 0 时 ， 播 值 曲线 能 很 好 地 逼近 于 被 插 函 数 . 在 实际 应 用 过 程 
中 ， 这 种 逼近 条 件 还 可 以 适当 地 放宽 ， 比 如 ， 插值 区 间 访 可 以 不 需要 很 小 ， 只 要 参数 取 值 得 
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当 ， 有 理 三 次 Hermite 播 值 样 条 也 能 很 好 地 逼近 被 插 函 数 .， 甚至 比 三 次 Hermite 插值 样 条 通 近 
效果 还 好 . 这 里 ， 我 们 先 给 出 “好 的 逼近 ”的 定义 . 
定义 3 设 [zwi,ziti] 上 的 三 次 Hermite 插 值 样 条 为 Hi(z), 有 理 三 次 Hermite 插 值 样 条 RH;(x)， 


RHe;= max |RHi(z)—-yw(x), 五 ci = 
1 


Ti<T<Ti+ 
则 当 RHe;: < Hei 时 ， Tk RH;(x) 比 Hi(x) 1E [zi zi+1] 上 对 y(x) AIREA. 
根据 定义 3， 可 由 不 等 式 RHe, < He, 解 得 在 每 个 区 间 上 参数 Xi 的 取 值 范围 ， 在 参数 X; 取 
值 范围 内 任 取 一 个 值 ,都 可 保证 在 整个 插值 区 间 上 生成 一 条 “好 的 逼近 ”曲线 . 
fll 设 被 插 函数 为 


max |Hi(z) — y(z)], 


Ti<T<Titl 


2 . /m 
y(x) = zg (22). 
节点 取 zi = ili = 1,2,… ,10)， 即 插值 区 间 h; = 1 (i= 1,2,… ,9)， 由 不 等 式 RHe; < Hei] 
计算 出 在 各 区 间 上 参数 X; 的 取 值 范围 表 1 给 出 了 对 于 给 定 的 插值 条 件 ，X; 的 取 值 与 各 区 间 的 
误差 最 大 值 . 图 2 在 同一 个 坐标 系 画 出 了 RH(x), H(z) 和 wy(z)， 实 际 上 三 者 在 区 间 [0,2] 上 几 
平 重合 ， 为 了 便于 区 别 观 察 ， 图 2 中 的 RH(z) 整体 上 升 0.04 个 单位 ， 互 (z) 整体 下 降 0.04 个 单 
位 .图 3 表示 RH (z) 5 H(z) 相对 于 y(x) 的 误差 曲线 ， 由 此 可 见 ， 只 要 参数 选择 的 得 当 ， 有 理 
三 次 样 条 插值 曲线 比 三 次 样 条 插值 曲线 相对 于 被 插值 曲线 有 好 的 逼近 ， 


表 1: ”对 给 定 的 插值 条 件 ， 和 i 的 取 值 与 各 区 间 的 误差 最 大 值 





Ti Vi di ^i RHeij Hei 
1 0.636619772 -0.636619772 1.742796037 0.000225101 0.001140174 
2 0.000000000 -0.500000000 1.855965883 0.000026309 0.002618395 
3 -0.212206590 0.070735530 1.980443075 0.000076469 0.000169431 
4 0.000000000 0.250000000 1.780004135 0.000062045 0.001940114 
5 0.127323954 -0.025464790 1.880158484 0.000004777 0.000548286 
6 0.000000000 -0.166666667 1.770215018 0.000046646 0.001325298 
7 -0.090945682 0.012992240 1.846344070 0.000007414 0.000554593 
8 0.000000000 0.125000000 1.769263948 0.000035726 0.000985197 
9 0.070735530 -0.007859503 1.829906858 0.000009867 0.000504587 
10 0.000000000 -0.100000000 
0 一 RH(X) 
esi 
图 2: 有 理 三 次 Hermite 插值 曲线 (上 )、 被 插曲 线 (中 ) 和 三 次 Hermite 插值 样 条 曲线 (下 ) 
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图 3: ”有理 三 次 Hermite 插值 曲线 (H) 和 三 次 Hermite 插值 曲线 ( 红 ) 的 误差 曲线 


采用 张 量 积 的 方法 ， 我 们 可 以 构建 分 片 有 理 三 次 Hermite 插值 曲面 . 有 理 三 次 Hermite 插值 
曲面 与 三 次 Hermite 插值 样 条 曲线 相似 的 性 质 . 

定义 4 VE f(z,y) Aog XE [a,b] x [cq 的 二 元 函数 ,，a = ro < x34 & < Em = b c= 
yo «yi € < yn 二 d 是 分 划 点 ， 记 


T— Ti " Y — Yj 


hi = ti+1 — Ti, hj = Yj+1 — Vj U= h ^ z h; 
i J 
称 
RFi(v) 
i:9i1 Hiato) 
RH(z, an = (RFi(u), RFi+1 (u), RG: (u), RGia(u))M RG(v) (8) 
RGir1(v) 


X [a,b] x [c,d 上 插值 于 f (m, y) 的 分 片 有 理 三 次 Hermite 插值 曲面 ， 其 中 


f(0,0)  f(01) fo(0,0) fo(0,1) 
f(1,0) FA,1) fl(1,0) fil1,1) 
fal0,0) fu(0,0) fu,(0,0) fu,(0,1) 
f40,0) fa) fall,0) full,1) 


给 定 函数 的 端点 、 端 点 的 一 阶 偏 导 数 及 二 阶 混合 偏 导数 的 值 ， 当 选择 适当 的 参数 时 ， 有 理 
三 次 Hermite 插 值 曲面 比 三 次 Hermite 插 值 曲面 有 好 的 允 近 . 


例 2 取 被 插 函 数 为 
f(z,y) = sin (22) COS (zv). 


分 划 点 为 4a=0<1<2=b,c= 一 1 <0<1=d， 即 插值 区 域 为 [0,2] x [-1,1,. 2428 f 
为 和 il = A2 = 1.7922， 可 计算 出 被 插 函 数 与 有 理 三 次 Hermite 插 值 曲面 在 [0,2] x [-1, 1] 上 的 最 
大 误差 为 3.8 x 10-“， 如 图 4. 而 被 插 函 数 与 三 次 Hermite 插 值 曲 面 在 [0,2] x [--1,1] 上 的 最 大 
误差 为 1.6 x 10-?， 如 图 5. 


第 3 期 谢 进 ， 等 : 有 理 三 次 Hermite 插值 样 条 及 其 逼近 性 质 391 

















i: 有 理 三 次 Hermite 播 值 曲面 、 三 次 Hermite 插值 曲面 及 被 插 函 数 的 图 象 几 乎 重合 ， 
为 了 观察 方便 ， 图 6 和 图 7 中 的 有 理 三 次 Hermite 播 值 曲面 与 三 次 Hermite 插值 曲面 分 别 下 降 
了 1 个 单位 . 





图 4: 有理 三 次 Hermite 曲面 与 被 插 函 数 的 误差 图 5: 三 次 Hermite 曲面 与 被 插 函 数 的 误差 





图 6: 有理 三 次 Hermite 曲面 (下) 与 被 插 函 数 (上 ) 图 7: Z% Hermite 曲面 (F) 与 被 插 函 数 (上 ) 


5 ”结论 


本 文 提 出 了 一 种 有 理 三 次 Hermite 插 值 曲线 . 利用 Peano-Kernel 定 理 ， 证 明了 其 逼近 性 
E. 只 要 参数 选择 得 当 ， 有 理 三 次 Hermite 插 值 曲线 比 标准 的 三 次 Hermite 插 值 曲线 具有 较 好 
的 逼 近 效 果 , 值得 一 提 的 是 ， 每 一 段 插值 曲线 只 含有 一 个 参数 ， 计 算 比 较 简 单 . 

如 果 选 择 参 数 的 合理 ， 这 种 有 理 三 次 Hermite 插值 曲 线 的 形状 具有 和 良好 的 可 约束 性 ， 可 另 
行文 研究 . 
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Rational Cubic Hermite Interpolating Spline 
and its Approximation Properties 
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2- Department of Mathematics and Physics, Hefei University, Hefei 230601; 
3- Department of Mathematics and Physics, Bengbu College, Bengbu 233000) 


Abstract: Improving the approximability of interpolating curves/surfaces is an important issue in 
the computer aided geometric design. A piecewise rational cubic Hermite interpolating spline with 
a single parameter is constructed in this paper. The approximation properties of the interpolating 
spline are studied and a method for improving the approximability of the curves/surfaces is intro- 
duced. Examples are given to illustrate the advantages of our method. The results show that, for 
given interpolating conditions, if the parameters are chosen properly, the introduced interpolating 
curves/surfaces can approximate the interpolated functions better than the standard cubic Hermite 
interpolating curves/surfaces. 

Keywords: rational cubic Hermite spline; cubic Hermite spline; Peano-Kernel theorem; shape param- 
eter; approximation 
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